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基于 DW 特征 参数 的 GNSS-R 海面 风速 反 演 研究 * 
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摘 要 : 传统 上 ， 机 载 风速 反 演 一 般 使 用 相关 功率 曲线 后 沿 直接 匹配 方法 ， 该 方法 需 从 实测 数据 中 提取 状态 信息 生成 
理论 波形 ， 而 从 实测 数据 中 提取 状态 信息 是 有 误差 的 ， 导 致 生成 的 理论 波形 存在 误差 ， 从 而 反 演 效果 不 好 。 对 此 ， 提 
出 一 种 新 的 基于 一 维 时 延 功率 波形 特征 参数 的 风速 反 演 方法 ， 该 方法 根据 风速 对 一 维 时 延 功 率 曲 线 时 延 窗 的 影响 ， 建 
立 以 一 维 时 延 功率 曲线 时 延 窗 作为 特征 参数 的 风速 反 演 模型 。 实 验 结果 表明 ， 新 方法 的 精度 为 0.77 ms， 相 比 传统 方 
法 ， 精 度 提 高 了 33.04%， 且 模型 简单 。 
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Research on sea wind speed inversion of GNSS-R based on DW feature parameters 
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Abstract: Traditionally, airborne sea surface wind speed inversion generally use relevant power curve directly matching method, 
this method need to extracted state information from the actual measurement data for generating theoretical waveform, however, 
the state information extracted from the actual measurement data is existed error, so that there is an error in the theoretical 
waveform, thus, the inversion result is not accurate. To this, this paper proposed a new method of wind speed inversion based 
on one-dimensional delay power waveform characteristic parameters, the method was based on the influence of wind speed on 
delay window of one dimensional delay power curve, it established the inversion model of wind speed by using the time delay 
window of one dimension delay power curve as the characteristic parameter. The experimental results show that the accuracy of 
the new method is 0.77 m/s. Compared to traditional methods, the accuracy is improved by 33.04%, and the model is simple. 
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0 ”引言 法 的 研究 工作 K3。 利 用 相关 功率 曲线 后 沿 直接 匹配 的 风速 反 演 
算法 是 传统 机 载 风速 反 演 的 主要 方法 ， 其 他 风 场 反 演 算法 是 在 
利用 GPS 散射 信号 进行 海面 高 度 测量 的 概念 自 1993 年 此 基础 上 发 展 起 来 的 外。 但 是 功率 曲线 后 沿 直接 匹配 算法 在 生 
Martion-Neria 博士 提出 以 来 ， 全 球 导航 卫星 系统 反射 信 


号 成 理论 模拟 波形 时 ， 需 从 实测 数据 中 提取 状态 信息 生成 理论 波 
逐 形 ， 而 从 实测 数据 中 提取 状态 信息 会 产生 误差 ， 这 将 导致 生成 

渐 发 展 成 一 个 新 的 分 支 上 0。 从 电磁 波 的 传播 理论 可 知 , 海面 反 的 理论 模拟 波形 存在 误差 ， 从 而 降低 了 反 演 精度 
射 信号 中 必然 含有 反射 面 的 相关 特征 信息 ， 其 中 包括 反射 信号 本 文 提出 一 种 新 的 基于 一 维 时 延 功 率 波形 (one- 
的 波形 和 极 化 特性 等 的 变化 。 因 此 ， 可 以 借助 对 反射 信号 的 接 dimensional delay power waveform，DW) 特征 参数 的 GNSS-R 
收 处 理 来 实现 反射 面 物理 特征 的 估计 并 进行 下 一 步 反 演 方面 的 海面 风速 反 演 方法 。 首 先 介 绍 了 GNSS-R 海面 风速 反 演 原理 ， 
相关 工作 器。 并 在 此 基础 上 对 传统 的 相关 功率 曲线 后 沿 直接 匹配 方法 进行 分 
GNSS-R 技术 在 测量 海面 风速 方面 取得 了 很 大 的 进展 ， 其 析 ; 然后 建立 起 以 一 维 时 延 功 率 曲线 时 延 窗 作为 特征 参数 的 海 
中 ， 机 载 方面 进行 的 飞行 实验 可 检验 性 强 ， 可 用 来 开展 反 演 方 面 风速 反 演 模型 ， 从 而 实现 了 海面 风速 高 精度 的 反 演 。 实 验 结 
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果 表 明 ， 与 传统 功率 曲线 后 沿 直 


新 反 演 方法 的 精度 提高 了 33.04%， 且 模型 简单 。 
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究 中 ，Zavortny 和 Voronovich 借助 于 双 基 
希 霍 夫 近 似 的 几何 光学 方法 提出 了 海 


GNSS-R 海面 风速 反 演 原理 


[ 接 匹 配 的 风速 反 演 方法 相 比 ， 


目前 ， 在 利用 GNSS 海面 散射 信号 进行 海面 风 场 反 演 的 研 


雷达 方程 ， 


勒 二 维 相 关 功 率 模型 中: 


G2(D) 


(YOF)=T Ho— x A 


4rR2(D)R2(D) 


看 散射 信号 的 时 延 一 多 普 


R+R, 


| 


C (1) 


x|S[f,(p)- fx6,(p)d’p 


这 是 一 个 二 重 积 分 ， 积 


菲 涅 耳 区 域 和 多 普 勒 区 的 交集 [9]。 其 中 : 
天 线 增益 ;Rk 和 分 别 表 示 GNSS 卫星 和 
距离 ， 和 表示 伪 码 相关 函数 ; 
在 给 定 的 卫星 仰角 和 接收 机 高 度 下 ， 随 着 风速 的 增加 ， 
直接 带动 海面 粗糙 度 的 增加 ， 准 镜 


区 域 为 散射 区 域 p ， 主 要 为 上 述 

: 7 为 积分 时 间 ; G 为 
接收 机 到 散射 区 域 的 
s 表示 多 普 勒 滤波 函数 。 


会 


古方 向 的 散射 信号 相关 功率 


减 小 , 从 而 互相 关 函 数 变 化 曲线 就 会 呈现 不 同 的 形状 四, 如 图 1 


所 示 。 
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由。 该 方法 首先 从 实测 


形 ， 


处 理 


态 信 


品 


如 图 


伶 场 


图 1 不 同 海 况 下 散射 信号 的 相关 功率 曲线 


因此 ， 可 以 通过 研究 海面 粗糙 度 变化 对 散射 信号 相关 功率 
的 影响 来 实现 海面 风速 的 反 演 。 


传统 的 GNSS-R 海面 风速 反 演 方法 


传统 海面 风速 反 演 是 利 / 
数据 中 提取 系统 状态 信息 ， 
化 处 理 ， 然 后 提取 实测 数据 波形 ， 


并 进行 归 


相关 功率 曲线 后 沿 直接 匹配 方法 


生成 理论 波 


进行 降 噪 


和 归 一 化 处 理 ， 最 后 将 归 一 化 后 的 实测 波形 与 根据 实测 状 
息 所 计算 出 来 的 归 一 化 后 的 理论 波形 进行 匹配 得 到 风速 ， 


2 所 示 。 
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过 


hinaXiv 例 作 期 刊 ， 


提取 系统 状态 信 ) 


理论 模型 模拟 波形 


归 一 化 


海面 风 场 
风 场 反 演 图 


CUO le t 


图 3 实验 场景 图 
利用 相关 功率 曲线 后 沿 直 接 匹 配方 法 反 演 风 速 。 反 演 结果 
如 图 4 所 示 。 其 中 红色 三 角形 线 为 实测 数据 ， 最 下 方 的 黄色 星 
线 为 风速 7.8 m/s 的 理论 曲线 ， 中 间 的 蓝 色 星 线 为 风速 8.8 m/s 
的 理论 曲线 , 最 上 方 的 绿色 星 线 为 风速 9.8 m/s 的 理论 曲线 。 通 
过 比较 可 见 ， 实 测 数据 与 8.8 m/s 的 理论 曲线 匹配 得 较 好 。 
GNSS-R 一 维 延 迟 图 像 ( 单 帧 》 
1 " J l 一 ”理论 78m/s 
加 二 
0.8 上 一 2 一 反 演 ‘8.8m/s 
07 上 
册 06 
加 05 
让 ol 
03 上 
0.2 上 
0.1 
0 4 1 
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4 ”风速 匹配 图 


最 后 ， 通 过 浮标 数据 对 风速 反 演 的 精度 进行 评估 。 图 5 所 


示 即 为 利用 


GNSS-R 反 演 得 到 的 风速 值 与 浮标 数据 的 对 比 结果 。 


2017 年 在 山东 威海 先后 进行 了 多 次 GNSS-R 机 载 飞 行 实 
验 ， 这 里 仅 对 某 一 次 机 载 飞行 实验 下 采集 的 数据 进行 处 理 。 实 


景 如 图 


3 所 示 。 


图 中 浮标 数据 显示 的 平均 风速 为 8.8 m/s; GNSS-R 实际 反 演 风 
速 为 7.4~11.5 m/s、 平 均 反 演 风 速 为 8.3 m/s， 由 此 得 到 其 风速 


反 演 精度 为 1.15 m/s。 
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第 五 架次 风速 反 演 结 果 ,RMSE=1.15m/s 
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图 5 与 浮标 数据 的 对 比 结果 


3 ”一 种 改进 的 基于 DW 特征 参数 的 GNSS-R 海面 

风速 反 演 方 法 

传统 机 载 风速 反 演 使 用 的 相关 功率 曲线 后 沿 直接 匹配 方法 
主要 存在 以 下 不 足 : 

根据 实测 数据 所 提取 的 状态 信息 生成 理论 波形 时 ， 由 于 从 
实测 数据 中 提取 状态 信息 会 产生 误差 ， 这 将 导致 在 生成 理论 波 
形 过 程 中 这 些 误 差 会 被 二 次 引入 ， 从 而 使 得 反 演 效果 不 好 ， 降 
低 了 反 演 精度 ， 而 且 计 算 量 大 ， 反 演 过 程 繁琐 。 

因此 ， 本 文 研究 提出 了 一 种 基于 DW 特征 参数 的 GNSS-R 
海面 风速 反 演 新 方法 。 
3.1 基于 DW 特征 参数 的 GNSS-R 海面 风速 反 演 机 理 

对 山东 威海 多 次 机 载 飞 行 实验 数据 进行 处 理 。 散 射 信号 的 
日 一 化 时 延 相关 功率 曲线 如 图 6 所 示 。 


0 1 2 3 要 的 好 70 6 
一 维 时 延 /0. l*chip 


图 6 不 同 风 速 的 散射 信号 归 一 化 功率 


图 6 可见， 风速 的 增加 ， 散 射 信号 对 应 的 时 延 值 向 右 仿 
移 ， 散 射 信号 时 延 相关 功率 向 时 延 更 大 的 方向 扩散 0 (图 中 的 
t 为 对 应 风速 的 时 延 窗 的 长 度 ) 。 因 此 ， 可 以 通过 研究 相应 风 
速 和 一 维 时 延 功率 波形 特征 参数 〈 时 延 窗 长 度 ) 的 关系 ， 即 建 
并 风速 与 时 延 窗 长 度 的 函数 关系 来 实现 海面 风速 反 演 。 
3.2 ”时 延 一 维 相关 功率 时 延 窗 长 度 的 提取 
提取 一 维 时 延 功率 曲线 时 延 窗 长 度 + 时 ， 为 了 更 好 地 反映 
波形 拖 尾 的 变化 趋势 ， 将 截取 阔 值 定义 为 el。 
本 文 研究 选取 了 三 组 风速 进行 分 析 ， 分 别 是 3.9 m/s、4.3 
m/s、8.8 m/s。 现 提取 这 三 个 风速 下 的 一 维 时 延 功率 波形 的 时 延 
窗 长 度 。 


刘 原 华 ， 
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3 
二 
D> 


等 : 基于 DW 特征 参数 的 GNSS-R 海面 风速 反 ; 


风速 为 3.9 m/s 时 , 对 应 的 时 延 一 维 相关 功率 波形 时 延 窗 的 


大 小 变化 如 图 7 所 示 。 
1040 TH | 
1020 | | 
2 1010 上 - | | 
由 | | 下 
妃 | | | | 
上 起 1000H | hw ' 
酝 | | | | 
东 eol | | 
黑 | i 
所 980- | | 
be | | | 
a on | | 
960 | | | 
0 10 20 30 40 50 60 
时 刻 
图 7 时 延 一 维 相关 功率 波形 时 延 窗 的 大 小 
由 图 7 可 见 ， 时 延 窗 的 大 小 在 1 000 ns 上 下 波动 ， 波 动 范 


逢 为 960 ns~1 040 ns。 
风速 为 4.3 m/s 时 , 对 应 的 时 延 一 维 相关 功率 波形 时 延 窗 的 


大 小 变化 如 图 8 所 示 。 
1080 关 一 
1070 |- | | | 
时 106o|- | || 
及 1050- | | 
| | | 
el bt | | | | | 
起 -OF 一 河 1 Lr T bf 
才 | | | | 1 
划 1oao|- | i | 
紧 | | bh 
人 1020|- | 由 | 
四 | rm 机 
各 1010|- | | 
| | | | 
1000 上 | | | P| a 
0 10 20 30 40 50 60 
时 刻 
图 8 时 延 一 维 相 关 功 率 波形 时 延 窗 的 大 小 
由 图 8 可 见 ， 时 延 窗 的 大 小 在 1 040 ns 上 下 波动 ， 波 动 范 


为 1 000 ns~1 080 ns。 
风速 为 8.8 m/s 时 , 对 应 的 时 延 一 维 相 关 功 率 波形 时 延 窗 的 


大 小 变化 如 图 9 所 示 。 
1240 T 7 
1220 -| | 
四 1200 | | | 
衣 1180 ~ | | 
全 | 
已 | | | 
起 1160|- 一 | | 
革 | | 
最 1140 - | | 
蛙 | 1 
各 1120 - ! 村 
则 | 
各 1100|- | 
i 
时 刻 
图 9 时 延 一 维 相关 功率 波形 时 延 窗 的 大 小 
由 图 9 可 见 ， 时 延 窗 的 大 小 在 1 160 ns 上 下 波动 ， 波 动 范 


居 为 1 080 ns~1 240 ns。 
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结合 图 7~9 可 以 看 出 ， 随 着 风速 逐渐 变 大 ， 对 应 的 一 维 时 


延 相关 功率 波形 时 延 窗 长 度 明显 变 大 ， 风 速 与 一 维 时 延 功 率 
形 时 延 窗 长 度 呈 现 正比 例 关 系 。 并 且 风 速 发 生变 化 时 ， 时 延 
长 度 变 化 明显 易于 区 分 ,在 精度 上 也 满足 反 演技 术 要 求 。 因 此 


这 为 利用 时 延 一 维 相关 功率 时 延 窗 长 度 反 演 海面 风 场 提供 了 
据 依 据 。 


3.3” 反 演 模型 建立 和 精度 分 析 

通过 3.9 m/s、4.3 m/s、8.8 m/s 这 三 组 风速 对 应 的 时 延 窗 
度 ， 结 合 风 速 与 时 延 窗 长 度 的 正比 例 变化 趋势 和 新 反 演 方法 
反 演 机 理 ， 利 用 数学 建 模 理论 知识 ， 对 时 延 一 维 相关 功率 时 


波 


om 


了 


数 


长 
的 
延 


窗 的 大 小 + 和 风速 上 数据 进行 建 模 ， 经 建 模 发 现 寡 函 数 模 型 
能 反映 出 风速 和 时 延 窗 长 度 的 正比 例 变化 关系 ， 且 误差 小 ， 
度 优 于 其 他 经 验 模型 。 对 此 ， 本 文 仅 讨论 采用 响 函 数 模 型 进 
建 模 。 


最 
精 
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u=ax*Tr +b*rt+c (2) 
其 中 : a、b、c 为 待定 系数 。 
现 对 时 延 一 维 相关 功率 时 延 窗 的 大 小 和 风速 进行 拟 合 ， 如 


图 10 所 示 。 拟 合计 算出 a= 1.42*105，b= -0.00462，c =6.01。 
11 
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图 10 时 延 一 维 相关 功率 波形 时 延 袜 的 大 小 和 风速 拟 合 图 
最 后 , 通过 浮标 数据 对 风速 反 演 的 精度 进行 评估 。 图 11 所 
示 即 为 利用 GNSS-R 反 演 得 到 的 风速 值 与 浮标 数据 的 对 比 结果 。 
从 图 中 可 以 看 出 ， 反 演 风速 与 浮标 风速 观测 数据 具有 很 好 的 
致 性 ， 由 此 得 到 其 风速 反 演 精度 为 0.77 m/s。 
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反 演 风速 与 浮标 数据 的 对 比 结果 


图 11 
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方法 ， 通 过 对 山东 威海 海域 实验 数据 进行 处 理 和 分 析 ， 建 立 了 
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结束 语 


本 文 介绍 了 


种 利 


] DW 特征 参数 进行 海面 风速 反 演 的 新 


以 


维 时 延 功 率 曲线 时 延 窗 作 为 特征 参数 的 海面 风速 反 演 模型 ， 


从 Tf 
参数 进行 海面 风速 反 演 的 新 方法 能 够 有 效 地 反映 出 风速 的 变化 ， 
反 演 精度 为 0.77 m/s， 与 传统 相关 功率 
相 比 ， 精 度 提 高 了 33.04%。 因 此 ， 新 反 演 方法 具 
型 简单 等 优点 ， 


T 


i 实现 了 海面 风速 的 高 精度 反 演 。 结 果 显 示 ， 利 用 DW 特征 


线 后 沿 直接 匹配 方法 
有 精度 高 、 模 
大 大 提高 了 工作 效率 ， 简 化 了 反 演 流程 。 
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